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Abstract 

Q^. To achieve high focal spot intensities in heavy ion fusion, the ion beam must be compressed 



O 



longitudinally by factors of ten to one hundred before it is focused onto the target. The longitudinal 
compression is achieved by imposing an initial velocity profile tilt on the drifting beam. In this 
paper, the problem of longitudinal drift compression of intense charged particle beams is solved 
(— I , analytically for the two important cases corresponding to a cold beam, and a pressure-dominated 
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beam, using a one-dimensional warm-fluid model describing the longitudinal beam dynamics. 



I. INTRODUCTION 

High energy ion accelerators, transport systems and storage rings [l|, 12, ISl y, |5| have a 
wide range of apphcations ranging from basic research in high energy and nuclear physics, to 
applications such as heavy ion fusion, spallation neutron sources, and nuclear waste transmu- 
tation. Of particular importance at the high beam currents and charge densities of interest 
for heavy ion fusion are the effects of the intense self-fields produced by the beam space 
charge and current on determining detailed equilibrium, stability, and transport properties. 
In general, a complete description of collective processes in intense charged particle beams 
is provided by the nonlinear Vlasov- Maxwell equations |l| for the self-consistent evolution of 
the beam distribution function, fb{x, p, t), and the self-generated electric and magnetic fields. 
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the effects of finite geometry and intense self-fields often make it difficult to obtain de- 
tailed predictions of beam equilibrium, stability, and transport properties based on the 
Vlasov-Maxwell equations. To overcome this complexity, considerable theoretical progress 
has also been made in the develo pm ent and application of one-dimensional Vlasov-Maxwell 
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to describe the longitudinal beam dynamics for 
a long coasting beam, with applications ranging from plasma echo excitations, to the in- 
vestigation of coherent soliton structures, both compressional and rarefactive (hole-like). 
Such one-dimensional Vlasov descriptions rely on using a geometric-factor (^f-factor) model 
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8J, |85| to incorporate the average effects of the transverse beam geom- 



etry and the surrounding wall structure. Despite the many successful applications of such 
one-dimensional Vlasov models to describe the longitudinal dynamics of long costing beams, 
even a one-dimensional Vlasov description is often too complicated to be analyzed in detail, 
except in some simple limiting cases. The next level of description of the nonlinear beam 
dynamics is provided by the macroscopic fluid equations, which correspond to the flrst three 
momentum moments of the Vlasov equation, with some particular closure scheme which 
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The usual closure assumption for 



relates the higher moments to the flrst three [86|, 

coUisionless plasma is provided by the adiabatic equation of state {ds/dt = 0, where s is 

the entropy per the unit volume), which expresses the thermal pressure as a function of the 



density. In one dimension, this is given by pX~^ = const. |86L l87l l88j |. where A is the hne 
density of beam particles and p is the hne pressure. Such a one-dimensional fluid model, 
combined with the adiabatic equation of state and a g-factor model for the average electric 
field, fits in the class of one-dimensional fluid problems which can be solved exactly (8Sy 
using the formalism described in Sec. II. 

In presently envisioned configurations for heavy ion fusion, multiple, high-current, heavy 
ion beams are focused to a small spot size onto the target capsule. To achieve a high- 
intensity beam focused onto the target, the beams are first accelerated and then compressed 



longitudinally. One of the possible ways to compress the beam longitudina 
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drift compression scheme illustrated in Fig. Q 90l 
scheme consists of two parts. In the first stage, an initial tilt in the longitudinal velocity 
profile —Vf{x) is imposed on the long charge bunch with some particular line density profile 
Af{x). After time Tghape, the beam line density evolves to a profile Aj„(x), with velocity profile 
—Vin{x). At this time, an additional velocity tilt Vin{x) — Vcomp{x) is imposed on the beam, 
and the beam is left with line density Aj„(x) and velocity tilt —Vcomp{x). We refer to this stage 
as the beam shaping stage. This stage requires beam manipulation (imposing the velocity 
tilt) and is done when the charge bunch is very long. At this stage, the longitudinal pressure 
and electric field are negligible, and the beam dynamics is governed by free convection. 
During this stage the beam may or may not be compressed. The purpose of this stage is 
to shape the beam line density profile into a certain intermediate line density profile Aj„(x) 
with velocity profile —Vcomp{x), such that during the next stage, which we refer to as the 
drift compression stage, after further axial drift during the time interval Tcomp (see Fig. Hj), 
the beam is compressed longitudinally until the space-charge force or the internal thermal 



pressure stops the longitudinal compression of the charge bunch 
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99|. At the point of maximum compression, the velocity tilt profile is completely 



removed and the beam is left with desired line density Ao(x). The final focus magnets then 
focus the beam onto the target, and the beam heats and compresses the target fuel. Stated 
this way, the longitudinal drift compression problem in the beam frame is equivalent to the 
time-reversed problem of the beam expanding into vacuum with zero initial velocity profile, 
and the specified initial line density profile Ao(x) which is desired at the time of maximum 
compression before final focusing. In this paper we employ a one- dimensional warm-fluid 
model with adiabatic equation of state to study this problem analytically for arbitrary final 
(maximally compressed) line density profiles Ao(x). 



We consider here the two separate cases corresponding to a cold beam, and a pressure- 
dominated beam. In the case of a cold beam, the internal thermal pressure is negligible, 
and the dynamics of the drift compression is governed by the self-generated electric field. 
In the case of a pressure-dominated beam, the self-generated electric field is negligible, and 
the beam compresses under the influence of the thermal pressure and the initial velocity tilt 
—Vcomp{x). One of the compression scenarios considered for heavy ion fusion is neutralized 
drift compression, where the beam propagates through a charge-neutralizing background 
plasma as it compresses longitudinally. In such a scenario, the beam is compressed only 
against the internal pressure. For simplicity, the present analysis is carried out in the beam 
frame where the particle motions are nonrelativistic. The final results can be then Lorentz 
transformed back to the laboratory frame, moving with axial velocity — VJ, = —^bC relative 
to the average motion of the particles in the beam frame. 

This paper is organized as follows. In Sec. II, we briefly describe the one- dimensional 
warm-fluid model equations and the formalism used for solving them analytically. In Sees. Ill 
and IV, the general solution for the expansion problem (the inverse to the drift compression 
stage problem) is obtained analytically for the two cases corresponding to a cold beam, and a 
pressure-dominated beam. In Sec. V the general solution to drift compression stage obtained 
in Sees. Ill and IV is illustrated by several examples. In Sec. VI we briefly discuss the beam 
shaping stage. 

II. THEORETICAL MODEL 



In the present analysis, we employ a one-dimensional warm-fluid model 86|, |87|, |88|, |91| to 



describe the longitudinal nonlinear beam d yna mics with average electric field given by the 



g-factor model with etE^ = —elgdX/dx |79|, |80|, |8l|, |82, |83|, |8^ |85| . For example, for a space- 



charge-dominated beam with flat-top density profile in the transverse plane, g — 2 ln{rw/rb) 
J85[- Here, A(x, t) is the line density, e^ is the charge of a beam particle, r^ is the conducting 
wall radius, and rf, is the beam radius. Generally, the beam radius, and therefore the g- 
factor, are functions of the line density and the external transverse focusing, and can change 
during the beam compression. In most of the drift compression scenarios it is preferable to 
maintain the beam radius (and therefore the g-factor) constant during the beam compression 
by adjusting the external transverse focusing [20, IHII IH2] ■ Therefore, in the present analysis, 
we assume that the g-factor is a constant. 



The macroscopic fluid equations for the hue density X{x, t), the average l ong itudinal beam 



velocity v{x,t), and the longitudinal line pressure p{x,t) are given by 

dX d 
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dt + 9^^^") 
dv d 

dt dx 



0, 

elg dX 1 dp 



dw 

dx ' 



(1) 
(2) 



rrii, dx rUhX dx 

where p = {po/Xl)X^ for a triple-adiabatic equation of state. Here we have introduced the 

effective potential w deflned by 

2 A c„ A , , 

"^ = ^^a; + YaI' (') 

where c^ = elgXo/mi, and c^ = 3po/2mbXQ are constants with dimensions of (speed)^, nib is 
the mass of a beam particle, and Aq and po are constants with the dimensions of line density 
and line pressure, respectively. 

For future application in Sees. Ill- VI, in the remainder of this section we summarize 
well-estabhshed theoretical technique developed in fluid mechanics JS2| that can be used to 
solve the nonlinear fluid equations (Q) and (j21). By introducing the velocity potential (p, 
where v = d(f)/dx, we can rewrite Eq. as 



d(h v"^ 

^ + T + " = °- 



(4) 



The full differential of 6 then becomes 



— — c/x + ^:—dt = vdx — \- w] dt. 

dx dt \ 2 



Next, following Landau and Lifshitz 



89|, we introduce the Legendre transform 



(i{/) = d{xv) — xdv — d 



t - + u. 



+ tdl — + w 



Introducing x = — xf — t{w + v'^/2), Eq. ^ can be expressed as 

/n.2 \ 



dx = —xdv + td\ \- w 



tdw + {vt — x)dv. 



(5) 



(6) 



(7) 



It follows from Eq. ((7j) that x can be considered as a function of the new independent 
variables (fjw), and that 

. dx . dx .„. 

t = --, x-vt = ---. (8) 

dw dv 

Therefore, if the function x is known as a function of its arguments (w, w), then Eq. (jHJ gives 
{v,w) as implicit functions of {x,t). 



To obtain the equation for x we rewrite Eq. (JJ as 



d{t,x) ^ d{t,x) ^ d{t,x) ' ^ ' 

where 

d{a,b) _ da db da db , , 

d{x, y) dx dy dy dx 
Assuming that v is not a definite function oi w [v ^ v{w)], we multiply Eq. (jHl) by 
d{t,x)/d{w,v) and use the multiplication property for determinants. This gives 

0[W,V) 0[W,V) 0[W,V) 

The case when f is a definite function of w in some region of the (x, t) plane corresponds to 
a simple wave and will be considered later. Because A = X{w), Eq. (jllj) reduces to 

dX dx dX dt dt 

dw dv dw dv dw 

Substituting Eq. (jS)) into Eq. (jT^. we obtain the equation for x 

1 rfA /9x _ S^xA S^x ^ Q 
A dw \ dw dv"^ J dw"^ 

Note that Eq. (fT^ is a linear partial differential equation for the function x(f,w). By 
introducing the effective sound speed defined by c^ = Xdw/dX, we can rewrite Eq. (fTSjl as 

^-^ + c^^ = 0, (14) 

dw dv'^ dw'^ 

where c^ is to be regarded as a function of w. Equation (jT^ together with Eq. (jHJ can be 

used to obtain the solution to the system of equations (Q) and (j21) everywhere in the (x, t) 

plane except in the regions corresponding to simple wave solutions where v = v{c) 89]. In 

this case, the Jacobian A = d{v,w)/d{x,t) vanishes identically. In deriving Eq. (fTTjl . we 

divided Eq. Q by this Jacobian, and the solution for which A = is not recovered. Thus, 

a simple wave solution cannot be recovered from the general equation ()13p. 

If f is a function of A only, as in a simple wave, we can rewrite Eqs. (P) and ^ as 



Since 



dX/dt / dx\ dv/dt / dx' 

dX/dx \9t I .' dv/dx \ dt 



we obtain from Eqs. (fT3j) and (jT^ 



' dx\ 



v + 



dw 
dv 



(17) 
fl8) 



However, because v = v{X), it follows that {dx/dt)\ = {dx/dt)y, so that Xdv/dX = dw/dv 
{dw/dX){dX/dv) = {c^/X)dX/dv, and therefore 



v = ± I -dX. 



(19) 



Next, we combine Eqs. (HIj), (HH)) and ^ to give (dx/dt)^ = {dx/dt)x = v + X{dv/dX) 
V ± c(f ) . Integrating with respect to t then gives [89| 



X = t[v ± c(f )] + /(f). 



(20) 



where /(f) is an arbitrary function of the velocity v determined from the initial conditions, 
and c(f ) is given by Eq. (fTn|) . 

Equations (fT^ and (f^UI) give the general solution for the simple wave. The two signs 
in Eqs. ()19|) and ()20j) correspond to the direction of wave propagation: (+) is for a wave 
propagating in the positive x direction, and (— ) is for a wave propagating in the negative x 
direction. 

It is also convenient to solve Eqs. (0) and Q using the method of characteristics [89]. 
Multiplying Eq. (P) by c/X, and then adding and subtracting from Eq. Q, and making use 
of the relation dw/dx = {dw /dX){dX/dx) = {c'^/X){dX/dx), we obtain 



dv cdX I dv c dX^ 



0. 



We now introduce the new unknown functions 

J+ = v+ j -dX, 



J_ 



A 



-dX, 



(21) 



(22) 



which are called ^emann invariants. In terms of J+ and J_, the equations of motion take 
the simple form 






J-, 



0, 






j_ 



0. 



(23) 



The differential operators acting on J_|_ and J_ are the operators for differentiation along the 

curves C+ and C_ ( called characteristics) in the {x,t) plane given by the equations 

dx dx 

dt ' dt 
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C4 



V — c. 



(24) 



The values of v and c at every point of the (x, t) plane are given by the values of the Riemann 
invariants J+ and J_ which are transported to this point along the C+ and C_ characteristics 
from the region where the values of J+ and J_ (and therefore v and c) are known. Equations 
(I22|) and (j24j) are very convenient for numerical solution of the equations of motion and also 
for general analysis of the flow. 

The totality of available space-time generally consists of the regions where either a simple 
wave solution [Eq. (PUJ) ] or a general solution [solution to Eq. ^^^ together with Eq. (jH))] 
is applicable. The boundary between the simple wave solution and the general solution, 
like any boundary between two analytically different solutions, is a characteristic JBil]. In 
solving particular problems (see Sec. V), the value of the function x(t>, w) on this boundary 
characteristic must be determined. The matching condition at the boundary between the 
simple wave solution and the general solution is obtained by substituting Eq. (jHJ for x and 
t into the equation for the simple wave [Eq. ()20|)]. This gives 

^±c^ + f{v) = Q. (25) 

ov ow 

Moreover, in the simple wave solution (and therefore on the boundary characteristic), we 

obtain dw/dv = {dw/d\){d\/dv) = {c^/\){dX/dv) = ±c. Substituting into Eq. (^3)) then 

gives 

dx dwdx ^ ^. . dX , r( ^ n l'^a\ 



Equation (PH|) can be integrated to give 

X = -Jfiv)dv, (27) 

which determines the required boundary value of x- 

III. GENERAL SOLUTION 

Here we consider two separate cases. For the case of a cold beam {po = and ci = 0) 
it follows that w = c^(A/Ao) [Eq. 0], and c^{w) = Xdw/dX = w. In the opposite limit 
where we can neglect the electric field compared to the thermal pressure, it follows that 
w = (Cp/2)(A/Ao)^, and c^{w) = Xdw/dX = 2w. Introducing in place of w the variable 
c = y/nw, where n = 1,2, we can rewrite Eq. (J14p as 

^^ + i^_4^^ = 0, forn = l, (28) 

oc^ c oc ov^ 

Q2^n=2 Q2^n=2 



8 



0, for n = 2. (29) 



Equation ()29|) is an ordinary wave equation whose general solution is 

x'^='{v,c) = h{c + v) + Mc-v), (30) 



where /i and /2 are arbitrary functions. To find the general solution to Eq. ()28|) we first 
Fourier transform with respect to the v dependence. This gives 

n2 71=1 1 f)^n=l 

^^ + -^-^k\r' = o. (31) 

ac^ c oc 
Equation (jHTj) is Bessel's equation of order zero which has two Hankel functions, Hq {2kc) 
and Hq '{2kc) = [Hq '{2kc)]*, as independent solutions. Here ( *) re presents complex conju- 
gate. Using the integral representation of the Hankel function |lOO|, 

H^'\x) = -2z 1^°^ 0Ldt, (32) 



the general solution to Eq. ((2H1) can be expressed as 

Jkj^ dtA{k) J-^__ +JjkJ^ dtB{k) J-^__ , (33) 

where A{k) and B{k) are arbitrary functions. Finally, we can rewrite Eq. (jHH|l as 

/OO fjf 

-^==[/i(tc - v/2) + hitc + v/2)], (34) 

where /i and /2 are arbitrary functions such that the integrals in Eq. plj) converge. Equation 
(IM|) provides the general solution to Eq. (J2H)). 

In the regions pertaining to the simple wave solution, Eq. (jl9|) gives the general relation 
between velocity and the density or sound speed in the wave. In the two cases considered 
here, A/Aq = {cjCgY for a cold beam (n = 2), and A/Aq = c/cp for a pressure-dominated 
beam (ra = 1). For these two cases 

V = ±c + a, n = 1, (35) 

v = ±2c + a, n = 2, (36) 

where a is a constant. Equations (jH^ and (jHUj) together with Eq. (plj) give a simple wave 
solution for the two cases considered in this section. From Eq. (j221); the corresponding 
Riemann invariants can be expressed as 

Jl=^ = v + c, JT^=v-c, (n = l), (37) 

Jl=^ = v + 2c, JT^ = v-2c, (n = 2). (38) 
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IV. GENERAL SOLUTION OF THE INITIAL VALUE PROBLEM 

In this section we make use of Eqs. fl3U|) and (jMj) to solve the initial value problem for the 
case of beam expansion into vacuum. The initial conditions for this problem are zero flow 
velocity at every point, fo(a;, 0) = 0, and prescribed density profile A(x, 0) = Ao(a;), which 
expresses the initial line density as a function of x. At some later time t = tf, the density and 
velocity profiles will be given by the functions X{x,tf) and v{x,tf) which are the solutions 
to Eqs. (H} and (0). Since the equations of motion [Eqs. ((T)) and Q] are time-reversible, the 
flow described by A(x, t) = X{x, tf — t) and v{x, t) = —v{x, tf — t) are also solutions to these 
equations with initial conditions v{x,0) = —v{x,tf) and A(x, 0) = X{x,tf). At time t = tf 
this flow has zero velocity profile {v = 0) and the density profile is given by the initial profile 
for the expansion problem, i.e., A(x,t/) = Ao(a;). 

To solve the initial value problem we assume that the density profile Ao(a;), or equivalently 
the sound velocity profile Cq{x), decreases monotonically to zero at the beam boundary 
X = ±xo, is an even function of x, and is an invertable function for a; > everywhere where 
the density is non-zero. Therefore, we assume that at t = the inverted profile Xo(c) is 
known. The condition that Co(x) decreases monotonically to zero at the beam boundary 
means that no rarefaction wave is launched from the boundary into the beam as it expands. 
We will treat the case with discontinues in Cq{x) at the beam boundary in one of the examples 
in Sec. V. Since we are interested in the time-reversed problem of beam compression, we 
assume that multi-valued flow does not form as the beam expands. This is equivalent to 
considering only initial density profiles with first derivative decreasing continuously from the 
beam center to the beam edge. This guarantees that the portions of the beam with smaller 
density accelerate faster than the portions with lower density, and as a result, the flow is 
never multi- valued. We will treat the case with multi- valued flow in one of the examples in 
Sec. V. The region of flow in the {x, t) plane and its boundaries are illustrated in Fig. |21 It is 
obvious that the flow is symmetric under reflection, x —>■ —x, and therefore we need only to 
solve the equations in the region x > 0. In general, there are four regions of flow. (See Figs. 
121 and 121) Each is separated from the others by two characteristics, the C_ characteristic 
P on which v = nc, and the C_|_ characteristic Q on which v + nc = hcq. The boundary 
conditions are given by 

f(0,t) = 0, c(x, 0) = Co(x), t'(x, 0) = 0, for |x| < Xq, c[xb{t),t\ = 0, (39) 

where xq = Xh{t = 0) is the initial beam half- width, and Xbit) is the coordinate of the beam 
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edge. 

As evident from Fig. |21 the flow at every point in Region I is brought to this point by 
the characteristics originating from the x-axis at t = 0. Hence, the flow in this region is 
fully determined by the boundary conditions c{x, 0) = Cq{x) and v{x, 0) = 0. The flow at 
every point in the Region II, which is separated from Region I by the C+ characteristic Q 
originating from the origin in the (x, t) plane, is brought to this point by the characteristics 
originating from the a; = line where f (0, t) = 0, and from Region I. Therefore, the boundary 
condition for Region II is given by f (0, t) = and by the flow on the separating characteristic 
Q. The flow at every point in Region III, which is adjacent to the beam edge and separated 
from Region I by the C_ characteristic P originating from the beam edge at t = in the 
(x,t) plane, is brought to this point by the characteristics originating from the beam edge 
where c[xb{t),t] = 0, and from Region I. Therefore, the boundary condition for Region III 
is given by c[xb{t),t] = and by the flow on the separating characteristic P. The flow at 
every point in Region IV, which is separated from Region II by the P characteristic and from 
Region III by the Q characteristic, is brought to this point by the characteristics originating 
from Region II and Region III. Therefore the boundary condition for Region IV is given by 
the flow on the separating characteristics P and Q. 

The function xi^^^c) ^^^ Eq. (jH} provide the map of the flow region in the {x,t) plane 
illustrated in Fig. HI to the (f , c) plane (Fig. E)). The region is a triangle (0 < c < Cq) limited 
from above by the C+ characteristic Q which is a straight line in the (f , c) plane since on 
this characteristic J+ = v + nc = nco = const {n = 1, 2). This mapping is not one-to-one. In 
fact Regions I and II and Regions III and IV in the {x, t) plane map into the same regions 
in the {v, c) plane, which means that in the {v, c) plane there will be four functions, x^ : X^^ i 
X^^^ and x^^ ^ defined inside the area depicted in Fig. El which map the depicted (f , c) region 
back into Regions I, II, III and IV in the (x, t) plane, respectively, by means of Eq. (jSJ. 

Since at t = 0, f(x,0) = and c(0,x) = co(x), by making use of Eq. (jH)) we obtain the 
boundary conditions for x'^l'^? c) in the (f , c) plane in Region I, which can be expressed as 

t = = [ ^ J = 0, or equivalently, ( ^ ) =0, (40) 

Since f(0,t) = Oata; = 0, the boundary condition for x^^('^? c) is 

(dx''\ 



dv I ^ 
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0. (42) 



The second boundary condition for x^^(^) c) reflects continuity of the mapping in Eq. 



dv I , \ dv I , 



(43) 



'dx''\ fdx'\ ^44^ 



dc , \ dc , 



By making use use of c[xb(t), t] = 0, the deflnition c^ = Xdw/dX = nw and Eq. (jSI), we obtain 
the boundary condition for x^^^(f , c), 



dc I ^ 

/ c=0 



0. (45) 



The second boundary condition for x^^^i'^-, c) reflects continuity of the mapping in Eq. 

(dx^\ ^ / V\ 
\ dv \ dv 

\ / v=nc \ / v=nc 

(dx^\ ^ (dx[\ 
[ dc [dc 

\ / v=nc \ / v=nc 

Finally, the boundary conditions for x^^(^? c) reflects continuity of the mapping in Eq. 

(dx^ ^ (dx^\ 

\ dv j ^ \ dv j ^ 

\ / v+nc=nco \ / v+nc=nco 



(46) 
(47) 



(4J 






(50) 



\ dc I , \ dc I , 

\ / v+nc=nco \ / v+nc=nco 

and 

(dx^ ^ (dx^\ 

\ dv \ dv 

\ / v=nc \ / v=nc 

m <^) ■ 

\ / v=nc \ / v=nc 

Next, we consider separately the two cases corresponding to a cold beam, and a pressure- 
dominated beam. 

A. Pressure-dominated beam 

The general solution for the case of a pressure-dominated beam is given by Eq. ()30|). To 
satisfy the boundary condition in Eq. (j4(Jj) . we are required to choose 

X' = f(c-v)-fic + v). (52) 
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Substituting Eq. (|^ into Eq. (HU, we obtain f'{c) = xo(c)/2, and therefore 

I re 

/(c) = 7; Xo{c)dc. (53) 

Here, we have chosen the integration constant so that /(cq) = 0. Substituting Eq. (|3^ into 
Eq. (j3^ then gives 

x' = - r ''xo{c)dc. (54) 

Z Jc+v 



In Region II, to satisfy the boundary conditions in Eq. (j42|) . we are required to choose 

X^^ = g{.c-v)+g{c + v). (55) 



Y / rc—v rc+v \ 

X^^ = - f y xo{c)dc + J xo{c)dc] . (56) 



To satisfy the boundary conditions in Eqs. (P3|) and (jiij) we choose g{c) = f{c). Hence, the 
solution in Region II is given by 

]^ / rc—v rc+v 

/ xo{c)dc+ / 

'Co Jco 

It is readily shown that the solution in Region HI which satisfies all of the boundary condi- 
tions in Eqs. (j^ . (jl6|l and (jTfjl is given by 

X''' = -I (£" 2;o(c)cic + ^^'^ xo{c)dc^ , (57) 

and the solution in Region IV which satisfies all of the boundary conditions in Eqs. (j48p . 
P?|l . p]|l and dSH) is given by 



X^^ = - / Xo(c)dc. (58) 

Z Jv—c 

Finally, using Eq. (jH)) and the definition c^ = Xdw/dX = 2w, we obtain the solutions in 
Region I, 

X - vt = -[xo{c- v) + xo{c + v)], 

1 

t = yM^ ^ ^) ~ ^o{c + v)], (59) 



in Region II, 



in Region HI, 



X — vt = -[xo(c — v) — Xo(c + f )], 

1 



X - vt = -[xo{v - c) + Xo{c + v)], 

t = yM^ ~ c) - Xo{c + v)], (61) 
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and in Region IV, 

X -vt = -[xo{v - c) - Xo(c + v)], 

1 

t = —[xo{v-c) + Xo{c + v)]. (62) 

Equations (jHn|) and (fHUjl give the implicit solution describing the expansion of a pressure- 
dominated beam. We can also obtain the formulas for the asymptotic solution as t ^ oo or 
c — > 0. Indeed, for t — ^ oo or c — > the flow is almost entirely in Region IV. Using Eq. ()62|) . 
we obtain 

_ Idx^ _ xojv) x"{v)c 3 

c dc ~ c ^ 2 +^^^^' 

x-vt = ^ = -x'o{v)c + Oic^). (63) 

Finally, in the leading approximation, we can rewrite Eq. ()63p as 

X(xA) 1 /x\ / X a; p /„,x 

^'^- :Xo{-), v{x,t) = -, fort^oo. (64) 



Ag CqL \t/ t 

Evidently, the density profile given by Eq. ()64j) is correctly normalized. 

The same solution can be also obtained from a kinetic description. It's been shown in 
Ref. [8a| that Eqs. (jT]) and (J2I) [together with the adiabatic pressure relation p = Poi^/^o)^] 
are the two key moments of the kinetic Vlasov equation for a waterbag distribution function 
( / = const, in an enclosed area of phase space). Indeed if we denote the upper curve in 
Fig. m as f+(x, t) and the lower curve as f_(x,t), than by multiplying the Vlasov equation 
for f{x,v^,t) 

by 1 and by f^, integrating over Vx, and keeping in mind that / = const inside the region 
limited from above by v+{x,t) and from below by f_(x, t), we obtain 

— K(x, t) - v^{x, t)] + 2 ^K(^' ^f - ^-(^' ^f] = 0' (66) 

^^^[v^(x,tY-v.ix,tf] + ^-^[v^ix,tr~v.{x,tY] = 0. (67) 

Introducing the line density and flow velocity defined by 

A(x,t) = [,^(o,0)-._(0.0)] ^'^^'^''^"''^'^''^^' ^''^ 

vix,t) = -[v+{x,t) -v_{x,t)], (69) 
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where Aq is the density at x = at t = 0, we can rewrite the Eqs. ()66p and ()67p in faniihar 
form 

-A(x,t) + —[X{x,t)v{x,t)] = 0, (70) 



|[A(x,t).(x,t)] + |:[A(x,t).(x,t)^] + ^ ^"-^^°'°^~;-^°'°^J' |:A(x,tf 



^0 

Comparing with Eq. Q, we obtain c^ = Adty/rfA = ([t;+(0, 0)-i;_(0, 0)]V4Ao)A^ or c(x, t) = 
{[v+{0,0) -v^{0,0)]/2Xo)X{x,t) = {l/2)[v+{x,t) -v^{x,t)]. Therefore, v+{x,t) = c{x,t) + 
v{x,t) and v^{x,t) = v{x,t) — c{x,t). If the initial profiles are given by v{x,0) = and 
c(x, 0) = Co(x), then t>+(x, 0) = Co{x) and f_(x, 0) = — Co(x). Since Eq. (|H^ represents 
the free-streaming motion of the particles in phase space along straight-line trajectories, we 
readily obtain the expressions for t>_(x,t) and v+{x,t), 

f + (x, t) = Co [x — v^ (x, t)t] , (72) 

w_(x, t) = ±co[x — f_(x, t)t]. (73) 

Here, the — sign in Eq. (ffHjl holds for |x| < xq (xq is the coordinate of the beam edge at 
t = 0) and corresponds to Regions I and II in Fig.|2l and the + sign holds for Xb(t) > \x\ > Xq 
and corresponds to Regions III and IV in Fig. |21 (see Fig. HJ. Equations (f72j) and (fTSj) can 
be rewritten as equations for c(x, t) and v{x,t), 

c(x, t) = - {co[x -{v + c)t] ±co[x-{v- c)t]} , (74) 

v{x, t) = - {co[x - {v + c)t] ip co[x - (t; - c)t]} . (75) 

By adding and subtracting Eqs. ()74|) and (f73j) . and inverting the resulting equations, we 
obtain 

±xo(f + c) = X — (f + c)t, (76) 

xo(c — f) = X — (t> — c)t, (77) 

for Regions I and II, and 

±xo{v + c) = X - {v + c)t, (78) 

xo{v — c) = X — {v — c)t, (79) 

for Regions III and IV. The + sign in Eqs. (f7Bj) - (f7^ corresponds to Region I [Eqs. (fTBjl and 
(f77j) ] and Region III [Eqs. (f7S|) and (jTUI)]. and the — sign corresponds to Region II [Eqs. (f7B|) 
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and (jTTj)] and Region IV [Eqs. (fTSj) and (f7^] (see Figs. |21 and IH). The ± signs appear here 
because we have assumed an even initial profile co{—x) = co{x). 

Finally, by adding and subtracting, Eqs. (f7n |) - (f7m) take the form shown in Eqs. (j^UjI - ljU^ . 



B. Cold beam 



Here we use the general solution in Eq. ()34|) to solve the same initial value problem as 
discussed in the previous section, applied now to the case of a cold beam. To satisfy the 
boundary condition in Eq. (j4Uj) we are required to choose 

dt 



x\v,c) 



Substituting Eq. (IHUj) into Eq. (HU, we obtain 

1 (dx'\ 



[f{tc+v/2)-f{tc-v/2)]. 



(80) 



Xo[C) = - 



2c V dc 



v=0 



dz df{z) 



v^ 



-2 dz 



^1) 



Equation (jRT|) can be inverted by using the integral Abel transform in Appendix A. This 
gives 



/(^) 



2 f^o qxo{q)dq 



TT Jz i/g^ — z'^ 



Q{z < Co) 



(82) 



where Cq = Co(x = 0), Q{z < a) is the Heaviside step-function, and ^ > 0. Note from Eq. 
(j8(jp that in Regions I and II, where v < 2c, the argument of / in Eq. (j8(jp is positive, and we 
can use the form of / defined in Eq. (j82j) . For f > 2c (Regions III and IV), the argument of 
the function under the integral in the first term in general solution in Eq. ()34|) can become 
negative. Next, we show that the solution of the form in Eq. (jHUj) with f{z) continued to the 
regions where 2; < as f{z) = f{—z), or f{z) = f{\z\), will satisfy the boundary conditions 
in Eqs. fli^ - (jT7|) . Indeed, by expanding / in a Taylor series in Eq. (jHUI) for c ^ 0, we obtain 



(co+v/2)/c dt 



V^ 



--f{tc-v/2) 



fi-vm + ^n-vm 



In 



'co + v/2' 



+ 0{c' 



and therefore 



X 



III 



/(t;/2)ln 



'co- v/2' 



f{-v/2) In 



'co + v/2" 



O 



cMn(^ 



(83) 



(84) 



Differentiating with respect to c, and taking the limit c — i> in Eq. 



we obtain 



'dx 



iir 



[f{-v/2) - f{v/2)] 



dc 



+ 



c^O 



cln 



Co 



(85) 
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provided f = f{\z\). It readily follows that the continuity conditions in Eqs. ()46|) - ()47|) are 
also satisfied. Therefore, the solutions in Regions I and III are given by 

where 



..[f(tc + v/2)-f{tc-v/2)]. 



(86) 



11 J\i 



2 f^o qxo{q)dq 



0(|z|<co). 



^q'2 _ ^1 
To obtain the solutions in Regions II and IV, we note that the function 

df 

:[/(tc + t;/2) + /(te-t;/2)]. 



^7) 



X^^'^>,c) = - 



satisfies the condition in Eq. (021) • Also, if we choose / as in Eq. (|H7jl . the second term (and 
its first derivatives) in both equations (J8U|) and (j88p is zero on the dividing characteristic 
2c + f = 2co, and therefore all of the continuity conditions in Eqs. ()48 p - (jM]) are also satisfied. 
Equations ()86|) - (jHSj) together with Eq. (jH)) give the formal solution of the expansion problem 
for the case of a cold beam. Finally, substituting Eq. (jHTj) into Eqs. (jHUj) and (|HHj). changing 
the order of integration, and performing the integrations, we obtain 



2c /■l+''/2c 

X{v,c)= / dqy/qxo{cq)K 

X^\v,c) = / dq^/qxoicq)K 

TT J1-d/2c 



8c /■'=o/c 



dqqxQ^cq) 



{vl2cf-{q-lf 

Aq 
{vl2cf-{q-lf 

Aq 



(89) 



-.K 



-K Jl+v/2c J(^q + 1)2 _ (^y/2cy 



iv/2cy 



{q + l)^-{v/2cy 



(90) 



2c rl+v/2c 

X^^\v,c) = / dqy/qxoicq)K 

n Jv/2c-i 

dqqxo{cq) 



{v/2cf~{q-lf 



Aq 



Ac r^ac-i 

^ ^0 ^{v/2cf - (g - 1) 



-.K 



Aq 



2c rl+v/2c 

X^^iv,c) = / dqy/qxoicq)K 

n Jv/2c-i 



dqqxo{cq) 



{v/2cy - (g - ly 
vl2cf-{q-lf' 



(91) 



Ac r^/sc-i 

^ ^0 ^{y/2cY - (g - 1)^ 



5C 



co/c 



dqqxo{cq) 



n JI+V/2C Uq ^ I 




{v/2c) 



(92) 



Here, K is the complete elliptic integral of the first kind |lOO|. In Sec. V, we illustrate the 
application of these solutions with several examples. 

We can also obtain the formulas for the asymptotic solution as t -^ oo or c ^ 0. Indeed, 
for t — > oo or c — > 0, the fiow is almost entirely in Regions II and IV. Using Eqs. (j^ and 
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X 



we obtain 



-f,c) 



f{v/2) + -nv/2) 



In 



'cg-W2)^ 



(93) 



and therefore 



c2 4 \ c2/ 



2c (9c 

X — vt = 



9x"' "■ _ m2) , /|\ 



9f 2 v^^y 

Finally, in the leading approximation, we can rewrite Eq. ()94|) as 



(94) 



X(x,t) xo f X 
~ :9 



Aq Cgi \ 2,tc^ 

where f{z) = cox^g^z / cq) and 

, , 2 rK\-') 

9{z) 



I, v(x,t) = —, for t — >• oo. 



(95) 



TT Jo 



Ao(x) — z^fix. 



(96) 



Here Ao(x) = Ao(a;/xo)/Ao is the scaled initial line density profile. One can readily verify 
that the density profile given by Eq. ^^ is correctly normalized, and that g{—z) = g{z). 

V. EXAMPLES WITH DIFFERENT INITIAL DENSITY PROFILES 

In this section, we apply the formalism developed in Sec. II-IV to several examples with 
different initial density profiles. 



A. Parabolic density profile 

As a first example, we consider here the case of an initial parabolic density profile for 
A(x, 0) = Ao(a;) with 

M£) . [i _ (iL)1 e(|.| < X.). (97) 

Aq [ \Xo/ 

1. Cold beam 

For a cold beam, c^ = Xdw/dX = c^(A/Ao). Substituting Eq. ^7\ into Eq. (j^ and 
integrating, we obtain 



/(^) 



CqXo 



1- 



co 



Q{z < Co) 



(98) 



18 



Next we substitute Eq. (|^ into the integral 



r{a,b) 



°° dt 



1 Vt^TTT 



f{tc-v/2) 



CqXq [^ dt 



a + 6)2 A Vt^^T 



{h-t){a + t) 



and define 






/+(a,6) 



V623T(6 - 2a) + (2a6 - 1) ln(6 + VW^^] 



°° dt 



where we have introduced new variables 



f{tc + v/2) = r{h,a), 



2c ' 
In term of the new variables, it follows that 

x'''''(a,6)=J+(a,6)-J-(a,6) = 

CqXq 



2co + v 2co - f 

= — , a = 



2c 



(99) 



(100) 



(101) 



- 2b) Va^ - 1 -{b- 2a) Jy^ - 1 + {2ab - 1) In 



(a + 6)2 

/^.^^(a,6) = -J+(a,6)-/-(a,6) 

CqXq 



a + \/a?^^ 



6+ V52^T 



(a + 6)^ 



- 26)Va2 - 1 + (6 - 2a) J y^ - I + (2a6 - 1) In 



+ V^^^T 



(102) 



(103) 



By introducing the scaled variables x = x/xq, i = tco/xo and x = x/^oCo, we can rewrite 
Eq. (jHI) as 



t 



X 



{a + bf 



(a" - b^] 



' — b— 
da db 

— b— 
da db 



X, 



X- 



(a + b) 



d__d_ 
db da 



X- 



(104) 
(105) 



Finally, substituting Eqs. (fTTl^ and (fmH|l into Eqs. (fTTl^ and (fmK|l . and using Eqs. dMI), 
we obtain after some lengthy algebra the solution in Regions I and III, 

il-v) 



t 



4c2 



1 +t;)2 -c2 - 



4c2 



1 — f )2 — c2 H — In 



1 1 + V+ Jil + vf - c 



4 l-i; + ,/(l-t;)2-c2' 



'7,2 



X 



C^ + f - v^)Jil + v)^ -c^ + {& -V - V^^x [l - v)^ - c 



\2 _ 7,2 



and in Regions II and IV, 



t =^^/(l-^)^-c2 + 



2c2 



4c2 



(106) 
(107) 



1 + t;)2 - c2 + - In ■ 



1 l+v + J{l + v)^-(? 



1 — f — J{1 — f )2 — C 



7,2 



X 



{(?^V- 1)2)^(1 + Vf -C^-{C^-V- V^)^{1 - V)^ - C 

2i2 



(108) 
(109) 
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Here we have introduced v = v/2co and c = c/cq. Equations ()106p - ()109p can be easily 
inverted. The result is 



c 
Co 



V 



Co 



l+/^\ 



1 



X 
Xo 



2/ 



I + /2 



+ 



2co 



A / 2/ 
Ac ll + /^ 



l + f 






2/ 



I + /2 



(110) 
(111) 
(112) 
(113) 



where < / < 1 is the solution of the transcendental equation 

- 1 - f ^ 1 

t = J; In/. (114) 

The solutions in Eqs. pi()j) - ()114|l describe the familiar self-similar solution JSO, ISlI, |22] for 
a parabolic density profile and is plotted in Fig. El Using Eqs. (j^ and (jHEI), we obtain the 
asymptotic solution as t — i> 00 

2' 



A(x,t) _ _2 
An "" 



3:0 

2cnt 



1 



X 



2tCr 



v{x,t) 



X 

i' 



for t -^ 00. 



(115) 



The exact solution given by Eq. pi3|) and asymptotic solution given by Eq. ()115|) are com- 
pared in Fig. ini (line b) for CqI/xq = 50. 



2. Pressure- dominated beam 



For a pressure-dominated beam, c^ = Xdw/dX = Cp(A/Ao)^ and threrefore the initial 



profile for Xo(c) is given by xo(c)/a;o = Jl — c/cq. Using Eq. (jHH) and Eqs. (|7i|l and (f7H|l we 
obtain the implicit solution in Regions I and II, 



;;Nil2 



c + v = l — [x — {v + c)i\ 



7;\T\2 



c — f = l — [x — (f — c)t\ 



(116) 
(117) 



and in Regions III and IV, 



V — c=l — [x — {v — c)t\'^. 



(118) 
(119) 
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Solving Eqs. ()116|) - ()119|) for c{x,t) and v{x,t), we obtain 

1 



c = ^[Vl + 4t2 + 4xt + VI + 4P - m- 2], 



T + Ip [\/^ + 4^^ " 4^^ - \/l + 4t2 + 4xt\, 



(120) 
(121) 



for Regions I and II (0 < a; < 1, < t), and 

1 



c=^Vl + 4t2-4xt, 

_ X 1 
' = I-2P^ 



(122) 
(123) 



for Regions III and IV (1 < x, 1/2 < t). The solutions [Eqs. (|T^ - (|T^ ] are illustrated in 
Fig. [7| Using Eq. (jMj) we obtain the asymptotic solution as t -^ oo 

A(x,t) 



Ao 



c=^Jl -, v{x,t) = -, ioTt^oo. 

Cot V Cot t 



(124) 



The exact solution given by Eq. (jl22p and asymptotic solution given by Eq. (jl24|) are com- 
pared in Fig. IHl (line b) for CqI/xq = 10. 



B. Linear density profile 

The next example we consider corresponds to the initial linear density profile 

Ao(x) 



An 



1- 



X 
Xo 



Q{\x\ <xo). 



(125) 



1. Cold beam 



Here we repeat the intermediate steps in the previous example. For a cold beam, c^ 
Xdw/dX = c^(A/Ao). Substituting Eq. ()125p into Eq. (jUUj) and integrating, we obtain 

3/2 



m 



AcoXo 
3tt 



Co 



Q{z < Co). 



(126) 



Next we substitute Eq. 



into the integral 



I-ia,b)=l ^=|L=/(tc-^;/2) = |^-^^^^=|L=[(6-t)(a + t)p/^ (127) 



and define 



1 Vt^^I' 
I+{a,b) 



37r(a + 6)3 7i y^W^^ 



°° dt 



1 vw^^ 



fitc + v/2) = I-{b,a), 



(128) 
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where a and b are introduced in Eq. (llUlj) . By introducing the new integration variable a 
defined by 



sm a 



\{t + i){b-iy 



(129) 



the integral in Eq. p27|l can be expressed as 



/-(a,6)_64 /6^T[(6-l)(a + l)]3/2 f^/2 cos\a)[l - kY sm\a)f/^ 



XqCq 



SnV b+1 



(a + bf 



[l-Fsin'(a)]^ 



(130) 



where fc^ = (6 — l)/(6+l) andp^ 



■ l)/(a+ 1). The integral in Eq. (jlHOj) can be expressed 



in terms of elliptic integrals. Finally, using the definitions x 



I, III 



I {b, a) — I (a, b) and 



X 



II, IV 



'I (a, b) — I {b, a), and using Eqs. (jlOll) . p04|) . and ()105|) . we obtain after some 



lengthy algebra the solution in Regions I and III, 



vc 



■K 



1 + cY -t;2 



n 



l+v-c (1-c) 



+ nfl^^.(i^C^I-A' 



l + V + C (l + C^2_,-,2 
1-C 






(131) 



1 — f + c' (1 + c)^ — f ^ 



'1 + C)2 - V^ 



2v, 



A(l + 2v^-c^)c 



X 



TX 



n 



c il-cY-v^' 



+ n 



1 + c)2 -v"^ . 



1 — f — c (1 — c) 



V + c (1 + c 



n\2 — 1)2 



(132) 



.^2 _ ^2^ 



1 — t" + c' (1 + c)2 — f 2 



ir 



(l + c)2-t;2^ 



1 + 2t;^ - c' 



and in Regions II and IV, 

4 1 



3vrc2,/(l + c)2-t;2 



6c^w<^ n 



'l + t;-c (1 -cY -v"^^ 



7^\2 _ ^,2 



- n 



1 — V — c il — c) — v^ 



yl — f + c (l + c)2 — t>2yj 
((l + c)2-t;2)(2c2 + ^;2-l)E 
2 1 



1 +t; + c 1 + c 



2c(2c^ + v^ - ?>cv^ + c-l)K 



[I - cY - v^ 



(133) 



c) 



2_^2^ 



V(l + c)2-t;2 



'1 + C)2 - V^ 



X 



37^C2J(1 + C)2- 



f ^ 



-2 



6c'' 1 - c^ + 2r w n 



„-i2n 



■l + t;-c (l-c)2 -t;2 



- n 



1 — V — c (1 — c) 



^^2 .^2^ 



7^\2 „-.2 



\1 — f + C (1 + Cj^— f^ 

+ v{{l + cf - v^){4: + c^ - Av^)E 



l + t; + c' (l + c)2-v2y 

- 2ct;(4 - 4t;2 + c(2 + c - 3c^ + 6v^))K 
/(l-c)2-^;2^ - 



(134) 



;i-c) 



^^2_^2^ 



'1 + C)2 - V^ 



{1 + cy-v^j 

where v = v/2co, c = c/cq, x = x/xq and i = tco/xo- Here, K, E an d 11 are the complete 



elliptic integrals of the first, second, and the third kinds, respectively 100|. The solution to 
the expansion problem for the linear density profile Eq. ()125|) is given by Eqs. p31|) - ()134j) 
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and is illustrated in Fig. |H1 Using Eqs. ()131|) and ()132p . we can rewrite the solution in 
Regions I and III as 

(135) 

(136) 

This simple flow with uniform acceleration in Regions I and III exists only until t = tcr = 
2xq/cq. At this time the characteristic Q overtakes the edge of the beam, and the flow for 
t > tcr is entirely in Regions II and IV and is given by Eqs. (jl33p and (|134j) . Using Eqs. (jHSI) 
and (jl26j) we obtain the asymptotic solution as t ^ oo 

3/2 



V=^t, 




Xq 




\ ^2 


X 


Ao 2xo 


Xq 



\{x,t) 



4x 







37rcot 



X 
2tCn 



v{x,t) 



X 

i' 



for t -^ oo. 



(137) 



The exact solution given by Eqs. ()133p and ()134|) and asymptotic solution given by Eq. p37|) 
are compared in Fig. IHl (line a) for CqI/xq = 50. 



2. Pressure- dominated beam 



For a pressure-dominated beam, c^ = Xdw/dX = Cp(A/Ao)^. Hence, for the linear density 
profile in Eq. ()125l the initial profile for Xq{c) is 



Xq[c) 

Xq 



Co 



0(c < Co), 



(138) 



where xq is the initial beam half-length, and cq is the sound speed at x = 0. Substituting 
Eq. ()138|) into Eqs. (j3^ - (jU^ . we obtain the solution in Region I, 



1 — X — 1 — X — — 

c = ^, V = t ^, for < t < 1, t < X < I, 



1-t' 



l-P' 



(139) 



in Region II, 



X 



1 + t' 



l + t 



, for < t, < X < min(t, 1), 



(140) 



and in Region IV, 



t — X 



xt — 1 



, for 1 < t, 1 < X < t. 



(141) 



t2-l' t2_l 

Region III has disappeared. Indeed, using Eqs. (pT|) we obtain x = 1 and i = 1, so that 
the entire region consist of just one point. The corresponding solutions [Eqs. p39j) - ()141|l ] 
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are illustrated in Fig. QUI Using Eqs. ()64p and ()138|) we obtain the asymptotic solution as 

t ^ oo 

X(x,t) Xq / X \ , s. X , , , 

^'^ - — 1 , v{x,t) = -, for t^oo. (142) 



Ao Cot V Cot/ t 

The exact solution given by Eq. ()141|) and asymptotic solution given by Eq. ()142|) are com- 
pared in Fig. IHl (line a) for Cot/xo = 10. 

C. Flat-top density profile 

As discussed in Sec II, the general flow consists of regions where either a simple wave 
solution [Eq. (p?|) ] or a general solution [solution of Eq. (fT^ together with Eq. ©] is appli- 
cable. Up to now we have considered flows with no simple wave regions. This is guaranteed 
provided the initial density profile is smooth everywhere except at the beam edge where the 
sound speed is zero, co(xo) = 0. Here, we consider an example with an initial flat-top density 
profile, and a corresponding discontinuity in density at the beam edge, i.e., 

^ = e(|x|<xo). (143) 

This problem is equivalent to the one-dimensional expansion of a uniform-density gas into 
vacuum in a vessel in which the end walls are instantly removed. The (x, t) plane is shown in 
Fig. ^2 The flow consists of three regions. In Region I the gas is at rest, and the information 
that the walls have been removed, which is carried by the C_ characteristic P into the gas, 
does not reach this region. Region II is the region occupied by a simple rarefaction wave 
which is centered at t = and x = Xq in the [x, t) plane, and is described by the simple wave 
solution in Eq. ()2U|) . Region III is the region of interference of this wave and and its reflection 
from the origin (or another rarefaction wave coming from the other end of the gas region). 
This region is described by the general solution in Eq. (|T4J) together with Eq. (jH}. Regions 
II and III are separated by the C+ characteristic Q. On this characteristic the boundary 
condition in Eq. (f?fjl holds. To determine the function f{v) in Eqs. (f^n|) and (jTfjl . we note 
that for t = 0, X = Xq and therefore f{v) = Xq = const. Also note that the characteristics 
C+ bring the value of the Rieman invariant J^ = v + nc = nco = const to all points of 
Region II. Therefore, the solution in Region II is given by 

n 1 X — Xq /-, , .\ 

-Co ^^-^, (144) 



n+1 n+1 t 

n n X — Xq ,-,a^\ 

Co + , , " , (145) 



n+1 n+1 t 
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where n = 1,2, and the boundary condition for x on the separating characteristic Q where 
V + nc = nco is given by 

X\v+nc=nco = -XoV ■ (146) 

The second boundary condition [Eq. ()42|) ] is given by 

1. Cold beam 

The function x satisfying the boundary condition in Eq. p47|) is given by 

/OO fif 

-^==[/(tc + v/2) + f{tc - v/2)]. (148) 

Using Eqs. p46|) and ()148|) . we obtain the integral equation for the function /, 



OO 



dt 
xov = J^ -^==^{/[(co - v/2)t - v/2] + /[(Co - v/2)t + v/2]}. (149) 

By setting f = in Eq. ()149p we note that the function / has the form f{t) = g{t)Q{t < cq). 
Substituting this expression for f(t) into Eq. p49p . and changing the integration variable to 
X = t(l — v/2cq) — v/2cq, we obtain the integral equation for the function g, 

fi dxg{cox) 



I I ^^ = vxq- (150) 

ll-v/co ^[x + 1)(X - [1 - v/Cq]) 

Changing the integration variables in Eq. (J15(J|) according to y = \/x + 1, and introducing 
the new function p(?/) = g\cQ{2y^ — 1)], we obtain an integral equation of the Abel type, 



where a = Jl — v/2cq. Equation ()15H) is easily solved using the Abel transform described 
in Appendix A. We obtain 



p{y) = coxo-yjl - y^. (152) 

71 * 



Finally, since g{x) = pU 1/2 + x/2cq), we obtain 



2 /, /x^2 
Next we substitute Eq. ()153|) into the integral 



f{x) = coxo^jl - (^) Q{x < Co). (153) 
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and define 



I+{a,b) 



dt 



-.f{tc + v/2) = r{b,a), 



(155) 



'1 V^-r 

where a and b are introduced in Eqs. (jlOll) . By introducing the new integration variable a 
defined by 



sm a 



\{t + i){b-iy 



(156) 



the integral in Eq. p54|) can be rewritten as 



/-(a,fe) _ 8 b-l[{b-l)ia + l)Y/^ /-/2 cos^(a)[l-A;Vsin^(a)]^/^ 
xoco ~7rV&+l (a + b) Jo [1 - k^ sm\a)]^ ' ^ ' 

where k"^ = {b — l)/{b + 1) and p^ = {a — l)/(a + 1). The integral in Eq. p57|) can be 
expressed in terms of elhptic integrals. Finally, using the definition x = —lib, a) — I^{a, b), 
and making use of Eqs. p01|) . ()104|) and ()105|) . we obtain after some lengthy algebra the 
solution in Region III, 



t 



X 



7;\'^\ 



7rc\/(l + cy-v^ 



7rc\/(l + cy-v^ 



[(1 + c) 



^[(1 + c 



v']E 



cY - v^ 



[1 + C)2 - V^ 



2cK 



cY - v^ 



[1 + C)2 - V^ 



+n 



c? - v^ 



1 + c 



f ' (1 + cY — v^ 




- vc{2 + c)K 



(158 
'\-cf-v''' 
:i + c)2-t;2_ 

(159) 



where v = v/2co, c = c/cq, x = x/xq and i = tco/xo- Here, K, E an d 11 are the complete 
elliptic integrals of the first, second, and the third kinds, respectively |lOO[. The solution in 
Region II is given by Eqs. ()144D and ()145|) with n = 2, i.e., 

a; — 1' 






t 

X — 1 

i 



(160) 
(161) 



Using Eqs. ()160p and ()16H) we can determine the trajectory of the beam edge and the 
characteristics Q and P. At the beam edge, c = 0, and therefore from Eq. ()16U|) we obtain 
Xb(t) = xo + 2cot. On the characteristic P, f = 0, and from Eq. (J16H) we obtain xp{t) = 
Xq — Cot. On the characteristic Q, dx/dt = v + c = 4co/3 + (x — Xo)/3t. Integrating this 
equation, we obtain 



Xgit) = Xo + Cot 2 



{cot/xoyr\ 



(162) 
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The solutions given by Eqs. p58p - p61|) are illustrated in Fig. [T21 Using Eqs. (j^ and p53|) 
we obtain the asymptotic solution in Region III [0 < x < XQ{t)] as t —>■ oo 



— c = \l—[- — , v[x,t) = —, tor t ^ oo. (163) 



Ao iTCot\l V2tco/ ' ' t 

The asymptotic solution in Region 11 [xQ{t) < x < Xbit)] is still given by Eqs. ()16Up and 
(I16H) . The exact solution given by Eqs. ()158|) - p6H) and asymptotic solution given by Eqs. 
fjl63p and ()16U|) are compared in Fig. IHl (line c) for CqI/xq = 50. 

2. Pressure-dominated beam 

The function x satisfying the boundary condition in Eq. p47p is given by 

X{v,c) = f{c-v)+f{c + v). (164) 

Using the boundary condition on the characteristic Q, it follows that x = —XqV for t; + c = Cq, 
and we obtain 

-vxo = fico-2v) + fico). (165) 

Substituting t> = into Eq. p65j) . we find that /(cq) = and therefore /(cq — 2v) = —vxq, 
or f{x) = (xo/2)(x — Cq). Using Eq. ()164|) . we obtain 

X{v,c) = xoic-co), (166) 

and using Eq. (jHl), we obtain t = xq/c and x — vt = 0. Therefore, the solution in Region 111 
is given by 

1 X 

c=^,v = ^, for 1 < t, <x< t- 1. (167) 

The solution in Region 11 is given by Eqs. ()144|) and ()145|) with n = 1, i.e., 

l/x-l\ l/x-l\,-, - , , 

c=-[l-^y v = -i^l + ^y |t-l|<x<l + t. (168) 

In Region 1, the gas is undisturbed: c = 1 and v = 0, for < x < 1 — t and t < 1. The 
characteristic Q is given by the equation Xgit) = i — 1, and the beam edge is given by 
Xb{t) = 1 + 1 . The solutions given by Eqs. ()167p and ()168p are illustrated in Fig. UHl Since 
X in Eq. p66|l is a linear function of c and is independent of v , the asymptotic (t -^ oo) 
solution coincides with the exact solution in Eq. (|167|) . 
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D. Continuous density profile (no sharp edge boundary) 

Up to this point we have considered flows which do not form shocks and therefore are 
time-reversible. For such flows, the compression problem is equivalent to the time-reversed 
expansion problem. In this section, we consider an example of fluid flow which forms shocks. 
For simplicity, we consider here a beam without sharp edges in which Ao(a;) decreases to zero 
monotonically as x -^ ±00. In particular, we consider the initial density profile given by 

M^ = I . (169) 

Ao cosh^ (x/xq) 

Expanding flows with initially smooth profiles extending to x — i> ±00 such as in Eq. (|169|) 
are entirely in Regions I and II. Indeed Regions I and II are separated from Regions III and 
IV by the C_ characteristic P with v — nc = —ucq = {n = 1,2). However, for profiles such 
as Eq. ()169p . Cq{x) > for all |x| < 00, and therefore the entire region |x| < 00 maps into 
Regions I and II in the {v, c) plane (see Fig. E)). As a result, the solution for the flow is given 
by Eqs. (j3^ and (plj) for pressure-dominated beams, and by Eqs. (jHUj) and (|HH|) . with f{z) 
defined in Eq. (j82|) for cold beams. 

1. Cold beam 

For a cold beam, c^ = Xdw/dX = c^(A/Ao), and therefore c/cq = l/cosh(a;/xo). Substi- 
tuting Eq. ()169p into Eq. (jHSl) and integrating, we obtain 

fiz) = xoco (l - ^ j 0(z < Co). (170) 

Substituting Eq. ()17()j) into the integral for /~(a, b) and integrating, we obtain 



00 



ria,b)=l -^.^f(tc-v/2) = -^[V¥^-b\nib + VW^)], (171) 

and 

-^=L=/(te + t;/2) = /-(6,a), (172) 

where a and b are introduced in Eq. (|1(J1|) . Using the definitions x^ = I^{b,a) — I~{a,b) 
and x^^ = ~I^{b,a) — I^{a,b), and making use of Eqs. pOlj) . p()4j) and ()1()5|) . we obtain 
the solution in Region I, 



^ 2c2 



;i+^)2_c2_^(l_^)2_c2 1 

X = V- ^ h - In 

c^ 2 
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(173) 


'l+v+^{l+v)^-c^' 


, (174) 


_i-v-J{i-vy-c'^_ 



and in Region II, 



2c2 



2-c2 



[l + vf-c' + Jil-vf-c^ 1 

X = V— ^:;r-^ 1 — In 

c^ 2 



l+vY-c 



y\2 _ ^2 



(175) 
(176) 



Equations ()173p - p76p can be partially inverted to give 

•jPX^ + X 
{PX^ - 1) 

where 



.^ . ?^!g|±lzl). (177) 



1 v^ 

X = 27^ ^ 27^ ^, (178) 

cosh [x — 2ft) sinh [x — 2vt) 

and A = A/Aq = c^. The solutions given by Eqs. ()177|) and ()178p are illustrated in Fig. IT^ 
As evident from Eq. ()177p . {dv/dX)t > in some regions, which means that the regions with 
higher density accelerated faster than the regions with lower density, and eventually multi- 
valued flow is formed (see Fig.ll4|l. This is unlike the previous examples where dv/dX)t < 
for all t, and there was no multi-valued flow. 

2. Pressure-dominated beam 

For a pressure-dominated beam, c^ = Xdw/dX = Cp(A/Ao)^. Using Eq. ()169|) and Eqs. 
(ffi|) and (f7S|) . we obtain the implicit solution in Regions I and II, 

2 [ cosh^[x — {v + c)t\ cosh^[x — {v — c)t] J 



1/ I L 

2 [ cosh^[x — {v + c)t\ cosh^[x — {v — c)t] 



^"(^'^) = ^ \ --.2r. .. , .^^ - „„..2r. .. ... } ■ (180) 



The solutions given by Eqs. ()179|) and ()180|1 are illustrated in Fig. Cni 

VI. BEAM SHAPING 

In this section we consider the beam shaping problem referred to in Sec. I. That is, given 
an initial line density profile Am(x) at time t = and final line density profile A/(x) at time 
t = Tshape, what are the initial and finial velocity profiles, Vi„(x) and Vf{x) respectively. 
The beam shaping stage is necessary to prepare the beam density profile for the final drift 
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compression discussed in previous sections. Here, as in previous sections, we analyze the 
time-reversed problem. Therefore, the initial density profile Aj„(x) for the time-reversed 
problem is given by Eqs. (|95|) and (j96|) for a cold beam, or by Eq. (I64|) for a pressure- 
dominated beam, and the final density profile Af{x) illustrated schematically in Fig. ^ 
During the beam shaping stage, the longitudinal pressure and electric field are negligible 
and the beam dynamics is governed by free convection decribed by 

dx dt dt nso^ 

dv dv dX 

Here, Eq. (jl82|l follows from multiplying Eq. |T^ by dw/d\, and is equivalent to Eq. ((T)). 

Equation p81|) implies that the function v is constant along the characteristic given by 

dx/dt = V, which therefore correspond to straight lines given by 

x = vt + f{v), or v{x,t) = V[x — v{x,t)t], (183) 

where v{x,0) = V{x) and V^[/(f)] = v. Equation p83|) gives a general solution to Eq. 
(118 1|) for the velocity profile v{x,t). Substituting Eq. ()183|1 into Eq. ()182|) . we obtain t = 
—/'(f) + q{v)/X, where q{v) is an arbitrary function of v. Using the initial condition that 
at t = 0, A = Xo{v), we obtain the solution to Eqs. |TT^ . A = Ao(t;)/[l + t/f'{v)]. Note that 
l/f'{v) = V'[f{v)] = V'{x - vt), and Ao(t;) = Ao[l^(x - vt)] = A{x - vt), where A(x) is the 
initial density profile as a function of x. Therefore, the solution to Eqs. ()181|) and ()182|) is 
given by [91 1 

v{x,t) = V[x-v{x,t)t], (184) 

_ A[x-v{x,t)t] 
^^''^^^-l+tV'ix-vix^m- ^^^^^ 

Setting t = Tshape and introducing new the function Uin{x) = x + TshapeVin{x), 
we can rewrite Eq. (I184J1 as Vf{x) = v{x,Tshape) = Vin[x - Vf{x)Tshape] = 
{Uin[x - Vf{x)T] -[x- Vf{x)Tshape\} /Tshape, OT cquivaleutly, X = Uin[x - Vf{x)Tshape]- Fi- 
nally, using the definition of the function Uin{x), we can rewrite Eqs. p84|l and p85|l in a 
compact and manifestly time-reversible form giving the finial formal solution to the beam 
shaping problem, i.e., 

Vfix) = ^^^^, VUx) = ^'f^~\ (186) 

-^ shape ^ shape 

Ui„(x) rx rx rUf{x) 

Af{u)du = / Ain{x)dx or / Af{x)dx = / Ain{u)du. (187) 

Jo Jo Jo 
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Here, Uf{x) = f/j„^(x) is the inverse of function Uin{x) such that Uf[Uin{x)] = x, and we 
have assumed that A/(— x) = A/(x) and Aj^{—x) = Am(x). Examples applying results in 
Eqs. dUni) and dHZI) can be found in Ref. |9l|. 

VII. CONCLUSIONS 

To summarize, we have studied the longitudinal drift compression of an intense charged 
particle beam using a one-dimensional warm-fluid model. We have reformulated the drift 
compression problem as the time-reversed expansion problem of the beam with arbitrary 
line density profile and zero velocity profile. We have obtained exact analytical solutions 
to the expansion problem for the two important cases corresponding to a cold beam, and 
a pressure-dominated beam, using a general formalism which reduces the system of warm- 
fluid equations to a linear second-order partial differential equation. We obtained simple 
approximate analytical formulas connecting the initial and final line density profile and flow 
velocity profile for these two cases. The asymptotic velocity profiles are linear in both cases, 
and correspond to free expansion as t —>■ cxd. The scaled density profile for a pressure- 
dominated beam far from the compression point was shown to be the functional inverse 
of the compressed density profile in Eq. (jMj) . For a cold beam, the profiles are connected 
by the Abel transform [Eqs. (jU3jl and (jHSI)]- The general solution has been illustrated for 
parabolic, linear, and fiat-top initial (compressed) line density profiles. For the case of a 
parabolic density profile, we have recovered the familiar self-similar solution [20, ISjJ, |22| • We 
have illustrated the formation of multi-valued flow with the exactly-solvable example in Eq. 
(I169|) . and identified the conditions for shock-free compression. 
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APPENDIX A: ABEL TRANSFORM 

Here we use the following definition of the Abel transform 

/(.) = A[giz)] = 2 r ^^^. (Al) 
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The inverse Abel transform is given by 



i-irx/^M 1 rdf{z) dz 



9(^) = ^-V(^)]---i ^^TTW- (^^' 



The fact that Eq. ()A2|) is indeed the inverse of the Abel transform in Eq. ()A1|) can be checked 
by direct substitution of Eq. ()A2|) into Eq. ()A1|) . Changing the order of integration leads to 

1 f°° 2xdx f°° df{t) dt 



Vr Jz y/x"^ - ^2 Jx dt y/t"^ - x"^ 

1 f°° , df(t) /•* 2xdx , , , 

Tl h dt Jz Vt^ - XWX"^ - Z^ 

Making the change of variables, sin(g) = \/x^ — z^/y/t"^ — z'^, the final integral in Eq. ()A3|) 
is evaluated to be 

TT, (A4) 



iz i/t^ _ x'^y/x'^ — z'^ 

and Eq. ()A3|1 becomes 

/(z) = - ^~ rft^ = f{z) - /(oo). (A5) 

Therefore, for functions /(z) such that fipo) = 0, Eq. ()A5j) is an identity. 
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FIGURE CAPTIONS 

Fig.l: Schematic of the two stages of drift compression described in Sec. I, the beam 
shaping stage and the drift compression stage. 

Fig. 2: The area in the {x, t) plane occupied by the four regions of flow. 

Fig. 3: The area in the (f , c) plane occupied by the four regions of flow. 

Fig. 4: Schematic in {x^v^) phase space of the flow for a pressure-dominated beam. 

Fig. 5: Plots of the normalized density A(a;,t)/Ao (a), and the flow velocity v/2cq 
(b), at CqI/xq = 1 as functions of x for a cold beam. The initial density profile (dotted line) 
is given by Eq. (jOTjl . 

Fig. 6: Plots of the normalized density A(x,t)/Ao as a function of x for a cold beam 
at CqI/xq = 50. The solid lines are the exact solutions. The dotted lines are the approximate 
solutions given by Eq. (jHH). The initial profiles are given by: (a) Eq. p25|) . (b) Eq. (P7|) . 
and (c) Eq. 



Fig. 7: Plots of the normalized density \{x, t)/Xo for cqI/xq = 0.5 (a) and for Cot/xo = 2 (c), 
and the flow velocity v/2co for cot/xo = 0.5 (b) and for cot/xo = 2 (d), as functions of x for 
a cold beam. The initial density profile (dotted line) is given by Eq. ()125p . 

Fig. 8: Plots of the normalized density A(a;,t)/Ao as a function of x for a pressure- 
dominated beam at cqI/xq = 10. The solid lines are the exact solutions. The dotted lines 
are the approximate solutions given by Eq. (jMj) . The initial profiles are given by: (a) Eq. 
dnSI), (b) Eq. (jnZ|), and (c) Eq. (fTIH|l . 

Fig. 9: Plots of the normalized density A(x,t)/Ao for Cot/xQ = 0.5 (a) and for Cot/xQ = 2 
(c), and the flow velocity v/cq for cot/xo = 0.5 (b) and for cqI/xq = 2 (d), as functions of 
X for a pressure-dominated beam. The initial density profile (dotted line) is given by Eq. (j97|) . 
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Fig. 10: Plots of the normalized density A(x, t)/Ao for Cot/xo = 0.5 (a) and for cot/xo = 2 
(c), and the flow velocity v/cq for cot/xo = 0.5 (b) and for cqI/xq = 2 (d), as functions of x 
for a pressure-dominated beam. The initial density profile (dotted line) is given by Eq. (jl25p . 

Fig. 11: The area in the (x,t) plane occupied by the three regions of flow. 

Fig. 12: Plots of the normalized density A(a;,t)/Ao for Cot/xo = 0.5 (a) and for CqI/xo = 3 
(c), and the flow velocity v/2co for Cot/xo = 0.5 (b) and for Cot/xo = 3 (d), as functions of 
X for a cold beam. The initial density profile (dotted line) is given by Eq. ()143|1 . 

Fig. 13: Plots of the normalized density A(x,t)/Ao for CqI/xq = 0.5 (a) and for Cot/xo = 2 
(c), and the flow velocity v/cq for cqI/xq = 0.5 (b) and for cqI/xq = 2 (d), as functions of x 
for a pressure-dominated beam. The initial density profile (dotted line) is given by Eq. 



Fig. 14: Plots of the normalized density A(a;,t)/Ao (a) and flow velocity f/2co (b) at 
cot/xo = 1.5 as functions of x for a cold beam. The initial density profile (dotted line) is 
given by Eq. ()169p . 

Fig. 15: Plots of the normalized density A(x, t)/Ao (a) and flow velocity v/cq (b) at 
cot/xo = 1.3 as functions of x for a pressure-dominated beam. The initial density profile 
(dotted line) is given by Eq. ()lf)9|l . 
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